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Abstract. We present results of simulations performed with the Geant4 software 
code of the effects of Galactic Cosmic Ray impacts on the photoconductor arrays of 
the PACS instrument. This instrument is part of the ESA-Iferschel payload, which 
will be launched in late 2007 and will operate at the Lagrangian L2 point of the Sun- 
Earth system. Both the Satellite plus the cryostat (the shield) and the detector act 
as source of secondary events, affecting the detector performance. Secondary event 
rates originated within the detector and from the shield are of comparable intensity. 
The impacts deposit energy on each photoconductor pixel but do not affect the 
behaviour of nearby pixels. These latter are hit with a probability always lower than 
7%. 

The energy deposited produces a spike which can be hundreds times larger than 
the noise. We then compare our simulations with proton irradiation tests carried 
out for one of the detector modules and follow the detector behaviour under 'real' 
conditions. 

Keywords: Instrumentation: detectors, photoconductor - Galaxy: cosmic rays - 
ISM: cosmic rays 



1. Introduction 

The European Space Agency Herschel satellite, scheduled for launch on 
August 2007 with an Ariane-5 rocket, will operate at the Lagrangian L2 
point of the Sun-Earth system (see ESA web page: www.rssd.esa.int/SA- 
general/Projects/Herschel). Herschel is the fourth cornerstone mission 
of ESA Horizon 2000 program and will perform imaging photometry 
and spectroscopy in the far infrared and submillimetre part of the 
spectrum. The Herschel payload consists of two cameras/medium res- 




© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



gegapaper2_rev.tex; 5/02/2008; 12:01; p.l 



2 



Bongardo et al. 



olution spectrometers (PACS and SPIRE) and a very high resolution 
heterodyne spectrometer (HIFI). 

The PACS (Photo-conductor Array Camera and Spectrometer) in- 
strument will perform efficient imaging and photometry in three wave- 
length bands within 60-300/Ltm and spectroscopy and spectroscopic map- 
ping with spectral resolution between 1000 and 2000 over the entire 
wavelength range (60-210/xm) or short segments. 
PACS is made of four sets of detectors: two Ge:Ga photoconductor 
arrays for the spectrometer part and two Si-bolometer arrays for the 
photometer part. On each instrument side each detector covers roughly 
half of the PACS bandwidth (see Poglitsch et al., 2004). More about 
Herschel - PACS can be found in the PACS Web pages: pacs.mpe- 
garching.mpg.de and pacs.ster.kuleuven.ac.be. 

It is well known that cosmic rays may influence strongly the detector 
behaviour in space. The performances of detectors such as ISOCAM 
and ISOPHOT on board of the Infrared Space Observatory (ISO) 
were largely affected for time periods long enough to corrupt a large 
amount of data. Our goal, in the present paper, is to exploit our present 
knowledge on detectors and cosmic environment to understand how the 
detector performances change from the expectations. 

In this paper we focus our attention on the GeGa photoconductor 
arrays (hereafter PhC) of the PACS instrument, following its response 
under simulated cosmic rays irradiation. The detectors are similar to 
those on board of ISO and we expect therefore that they will be strongly 
affected by cosmic ray hitting. 

In addition we compare proton irradiation tests, performed on a sin- 
gle module of the PhC, with Geant4 simulations and draw conclusions 
about their performances. 

The paper is organized as follows: in Sections 2 through 5, the used 
simulation toolkit, the input detector design, the Galactic Cosmic Ray 
spectra and the Physics List are outlined. In Section 6 we report the 
results. Section 7 deals with the irradiation tests which were performed 
at UCL-CRC on one detector module and the comparison with our 
simulations. In the final Section 8 we summarize the results reported 
in this draft. 



2. Geant4 Monte Carlo code 

Monte Carlo simulations were accomplished with the Geant4^ (here- 
after G4) toolkit, which simulates the passage of particles through 

^ Every simulation was run with the 7.0 patch 01 version, with CLHEP 1.8.1.0 
with a gcc 3.2.3 compiler 
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matter. G4 provides a complete set of tools for all the domains of de- 
tector simulation: Geometry, Tracking, Detector Response, Run, Event 
and Track management. Visualization and User Interface. An abundant 
set of Physics Processes handle the diverse interactions of particles 
with matter across a wide energy range, as required by the G4 multi- 
disciplinary nature; for many physical processes a choice of different 
models is available. For any further information see the Geant4 Home 
page at: geant4-web.cern.ch/geant4/- 

3. The Detector Design 

3.1. PACS AS A SPECTROMETER 

The PACS employs two PhC arrays (stressed/unstressed) to perform 
imaging line spectroscopy and imaging photometry in the 60 - 300 
wavelength band. In spectroscopy mode, it will image a field of about 
50 X 50 arcsec, resolved into 5x5 pixels, with an instantaneous spectral 
coverage of ~1500 km/s and a spectral resolution of ~175 kni/s with 
an expected sensitivity {5a in Ih) of 3 mJy or 2.5x10^^^ W/m^. 

The high stressed module is cooled down to 1.8 K (whereas the low 
stressed is cooled to 2.5 K), the Front End Electronic (FEE) to 4 K. All 
main metal components are made of a light metal alloy called Erg-Al 
(AlZn5,5MgCu). 

3.2. The Geometry 

Figures 1 and 2 show the modeled detectors with Geant4. 

The most important part of the PhC is the cavity. Inside the cavity 
there is the active part of the PhC, the pixel. Each pixel is made 
of Ge:Ga (in the high stress configuration Ga doping is 10-*^^ cm~^, 
that is an atom of Ga for about 10^ atoms of Ge; this negligible Ga 
concentration is not further considered here). The pixel is a simple box 
of 1x1x1.5 mm. Each pixel is sustained by two boxes ('cubes') made 
of CuBe2 (a Cu, Be and Pb alloy). On the external parts of these two 
last boxes there are two different cylinders made of CuBe2 (pedestal 1 
and 2). A steel cylinder (the kugel) is placed on the right-hand side of 
pedestal 1. A thinner and smaller cylinder is located in the left-hand 
side. The sequence of a single pixel is: kugel, pedestal 1, cube, pixel, 
cube and pedestal 2 with its thinner cylinder. Between two sequences 
of pixels there is an AI2O3 (Aluminum oxide) insulator. The sequence 
of the 16 pixels starts with an insulator and ends with a kugel. The 
pixels have a single face that is not shadowed by the cavity, that is, 
clearly, the one facing the fore optics (see below). 
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The fore optics (hereafter FO) is the most difficult part to be de- 
signed with G4. Actually it is constructed as a unique block made 
of Erg-Al, from which they subtracted a cone non orthogonal to the 
respective pixel. We designed each cone individually, making use of 
the boolean solids. All the FO should be gold coated with a 8.95 g 
gold mass. Simple calculations of the surface area of the FO leads to a 
coating thickness of p»1.6 /xm. We did not reproduce such a coating. 

The FEE is made of two different volumes: a box and a trapezoid, 
both 0.5 mm thick. The Erg-Al components of the module have three 
different thicknesses. The small cubes sustaining the FO are 1.75 mm 
thick; the end of the horseshoe is 2.3 mm thick; the horseshoe is 4.1 
mm thick. 

Figure 1 shows the geometry of the PACS PhC (the single stressed 
and unstressed modules), as drawn in G4. Figure 2 shows the whole 
stressed PhC detector. This geometry is that used by the G4 software 
for the simulations. 

For the Galactic Cosmic Rays simulations we considered all the 
components 'external' to the detector, i.e. the PACS box, the cryostat 
and the spacecraft and called the ensemble 'shield' and approximated 
it with an Aluminum sphere, with thickness of 11 mm and at a temper- 
ature of 80 K, which is the nominal temperature of the passive cooled 
Herschel telescope. The inner radius is 20 cm and the outer 21.1 cm. 

For the irradiation tests the external environment was modelled as 
close as possible to the real one and it is described in details in §7. 

The G4 software then requires the design of a 'world volume' con- 
taining all the detector part and the specification of the environment 
in which the detector is placed. We made the world volume a little bit 
larger than the shield, i.e. is a sphere with a radius 21.2 cm. 



4. Galactic Cosmic Rays spectra 

Herschel will be injected in a orbit around L2. The L2 environment, 
as regards the Galactic Cosmic Ray (hereafter GCR) rates, has been 
measured by the WMAP satellite (Prantzos 2004) . This environment is 
similar to that of the geostationary orbit and therefore the GCR spectra 
are known. In the limit of the distance Earth-L2 the GCRs flux gradient 
is negligible. Therefore we use GCR fluxes of a geostationary orbit as a 
reasonable approximation. In particular we use the Cosmic Ray Effects 
on Micro-Electronics (CREME) model set both for the solar-minimum 
(GCR maximum) and for the solar-maximum environment (Nymmik et 
al. 1992). This model uses numerical models of the ionizing radiation 
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Figure 1, The PACS high (loft) and low (right) stressed module (left). 
Magenta: the fore-optics. Yellow: the gold coating of the FEE and the active pixel. 
Red: the CuBe2 contacts. Blue: the AI2O3 insulator. Gray: the cylindrical steel 
segment. Cyan: Erg-Al components. 

environment in near-Earth orbits to evaluate the resulting radiation 
effects on electronic systems in space^. 

A very diffuse background, called galactic vacuum by G4, which is 
a statistical view of the real interactions and diffusions of the Galactic 
Cosmic rays sources with and from the interstellar material, is defined 
and used as a good approximation of an empty space where an isotropic 
cosmic ray flux is present. This environment has the following physical 
properties: 

A = 1.01 g/mole (atomic mass) 

Z = 1 (atomic number) 

p = lxlO~^^ g/cm^ (density) 

P = lxlO~^^ Pa (pressure) 

Our simulation considers the most common GCRs in the near orbit 
Earth hcliospherc: protons, alpha particles and nuclei of Li, C, N, O 

^ The CREME version used in this work is CREME94, which was more exten- 
sively tested and compared with measurements, however there is an updated version, 
CREME96, with a web page: https://creme96.nrl.navy.mil/ 
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Figure 2. The PACS array. It contains 25 modules of 16 pixels each, making an 
array of 25x16. 

(sec Figure 3). It cannot include heavier nuclei, such as iron, since Fe 
physics has not yet been implemented in the Ion Physics of G4 (see 
§5). In Table I we report the rates (number of particles per second) of 
each GCR type computed by integrating the fluxes reported in figure 
3 over the shield surface and particle energy. Herschel will reach the 
L2 point in 2007 at a solar minimum but it will operate until an epoch 
of increasing solar activity. Therefore, we consider both solar activity 
phases. 

4.1. The Generate Source Particle Module 

G4 is built in such a way that once the geometry (volumes) and the 
physics are chosen and setup, the tool may be played from outside 
the main program through an User Interface (UI) easy accessible and 
modifiable (i.e. a macro file). In this macro file the user specifies: (a) 
the graphical visualization and/or (b) the type and the number of 
particle hitting the detector and their energy. One single particle with 
a determined energy represents an event. 
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Table I. Rates (number of nuclei per sec- 
ond, on a sphere of 21.2 cm radius) for 
different GCR types and solar activity. 



GCR Type 


Solar min 
rates (#/s) 


Solar max 
rates (#/s) 


H 


11275.89 


4553.28 


He 


1102.56 


458.67 


Li 


7.46 


2.96 


C 


33.52 


13.92 


N 


8.96 


3.68 







13.00 




Figure 3. The differential GCR SEDs at the solar niiniiinini (left) and maxi- 
mum (right) . Eaeh curve eorresponds to a type of particle (from top to bottom: 
proton, a, C, O, Li, N). The rates reported in Table I are the integrals over 
energy and refer to the full shield surface. 



We simulated the GCR seedings via the General Source Particle 
Module (GSPM; Ferguson 2000a, b,c), a code for the spacecraft ra- 
diation shielding, that is based on the radiation transport code G4. 
The user must specify the input parameters of the particle source: (1) 
the spatial and (2) the energy distribution. Every time G4 requires a 
new event, the GSPM randomly generates it according to the specified 
distributions. We adopt an angular distribution that is isotropic over 
a spherical surface (the vacuum sphere), without any preferred parti- 
cle momentum. Wc then set the energy distribution of each particle 
according to the SEDs shown in Figure 3 (see also §5.2). 
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5. The Physics Model 

We had to define a list of physical interactions with the shielding 
and the different components of the photoconductors in order to best 
reproduce the impact of GCRs on the photoconductors. The list is 
comprehensive of all the possible physical interactions: 

— Electro-magnetic physics: photo electric and Compton effect; pair 
production; electron and positron multiple scattering, ionization, 
bremsstrahlung and synchrotron; positron annihilation. 

— General physics: decay processes. 

— Hadron physics: neutron and proton elastic, fission, capture and in- 
elastic processes; photon, electron and positron nuclear processes; 
7r+, 7r~, K+, K^, proton, anti-proton, a"'", a~ , anti-a"'", anti-a", 
X'^, x" 1 anti-x'*', anti-x~, w"*", a;~, anti-w''" and anti-w" multiple 
scattering and ionization. We added the High Precision Neutron 
datasct, that is valid for neutrons till 19.9 MeV. 

— Ion physics: multiple scattering, elastic process, ionization and 
low energy inelastic processes. We added inelastic scattering (J. P. 
Wellish, private communication) using the binary light ion reaction 
model and the Shen cross section (Shen et al. 1989). 

— Muon physics: jU''", pT multiple scattering, ionization, bremsstrahlung 
and pair production; ii~ capture at rest; r'^ and r~ multiple 
scattering and ionization. 

For 'ion physics' of aZp/ia-particles, protons, and H isotope nuclei 
we used the cross section by Tripathi et al. (2000), holding for energies 
of up to 20 GeV per nucleon. The cross sections used for the other 
ions (Shen et al. 1989) are valid only up to 10 GeV per nucleon. This 
means that we had to cut out high energy ions of the GCRs (whose flux 
is lower than 10~^ /m^/s/sr). The treatment of Ion Physics is strictly 
valid only for nuclei not heavier than C. However we have applied it 
also to simulate N and O, which should not behave much differently. 

The model called by our physics list is different depending on the en- 
ergy of the incoming particles. Our physics list is based on QGSP_BIC^ 
list (J. P. WcUish, private communication; see the Geant4 web page for 
details). Modifications to such a list are: 

1. We added the High Precision Neutron dataset, which is valid for 
neutrons up to 19.9 MeV. 

^ Quaxk Gluon String Physics_Binary Cascade 
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2. For neutrons and protons, we introduced^ the Bertini cascade (Bertini 
& Guthrie 1979) up to 90 MeV and made the BIG starting from 80 
MeV. 

3. QGSP is assumed vahd from 10, rather than 12, GeV. 

4. We have taken into account the gamma- and electro-nuclear reac- 
tions, so we considered the electromagnetic_GN physics list. 



5.1. Thresholds 

In order to avoid the infrared divergence, some electromagnetic pro- 
cesses require a threshold below which no secondary particles will be 
generated. Then, each particle has an associated production threshold, 
i.e. either a distance, or a range cut-off, which is converted to an energy 
for each material, which should be defined as an external require- 
ment. A process can produce the secondaries down to the recommended 
threshold, and by interrogating the geometry, or by realizing when the 
mass-to-energy conversion can occur, recognize when particles below 
the threshold have to be produced. 

We have set different production thresholds (hereafter cuts) for each 
geometrical region. 

We decided: 1) not to use very small cut values (below 1-2 //m), since 
this could affect the validity of physical models at such small steps 
(V. Ivantchenko, private communication); 2) to set cut values ~ 1/5 of 
the volume thickness. We identified groups of volumes with the same 
thicknesses: in particular we defined 6 different regions (see Table II). 
Results are robust since they do not depend on the cut values, both 
default and the chosen ones lead to similar outputs, with an increasing 
of the secondary event production. The results (energy deposition) 
among the different simulations were always comparable at la. On 
the one hand this was expected, since most of the volumes have still 
the same default cut value, 0.7 mm. On the other hand it is possible 
that Erg-Al thresholds have already achieved convergence at the default 
values, that is we are very slightly tuning them. Although results are 
pretty similar (sec Table III) those with our chosen cuts allow a better 
tracking inside each detector volume. 

* It is well known that (1) the BIC does not reproduce well data for energies 
below 50 MeV (Ivanchenko et al. 2003). In contrast, the Bertini cascade works well 
below 50 MeV for all but the lightest nuclei. And (2) the BIC for pions and kaons 
has not been tested yet (J. P. Wellish, private communication). 
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Table II. The five regions of the photo- 
conductor array and their specific cut 
values. 



Region 


Cut Value 


Thickness 


Shield 


2 mm 


11 mm 




Highjnm 


0.7 mm 


3.02 - 4.1 


mm 


Lowjnm 


0.3 mm 


0.73 - 2.1 


mm 


High_mic 


100 ^im 


300 - 600 




Low_mic_a 


5 /im 


20 /jm 




Low_mic_b 


3 fjxn 


10 /um 





Table III. Comparison between the average val- 
ues of the hit numbers in three distinct cases: 
theoretical computation (see §5.2), G4 simula- 
tion without cuts, G4 simulation with cuts. This 
computation was done for the high stressed PhC 
without shield and secondary events and for the 
pixel volume. 







Default Cuts 


Chosen Cuts 


Volume 


Th. Hit 


(^limpacts) 


(^^impacts) 


Pixel 


56.43 


50.45±10.34 


58.18±12.81 



5.2. Theoretical computation of GCR flux 

In order to disentangle the contributions from the different simulated 
parts of the instrument (detector and shield) we evaluated the impact of 
the GCR flux onto the 'naked' detector (without spacecraft shielding). 
In Table III we report the results from the simulated GCR flux hitting 
the detector only. 

We computed the number of impacts, nimpacts; as nimpacts = vr<I>5, 
where ^ is the incoming particle flux in cm~^ sr~^), and S is the 
total surface area of the volume considered. 

The hit average values in the two cases (with G4 default cut values 
and with our cuts) arc compared to those compute theoretically. 
Fixing specific cuts increases the error bars, but (nimpacts) comes closer 
to the expected mean value. The chosen cut values are a fine tuning of 
the Erg-Al cuts and, as a result, this leads to a very precise tracking 
inside the detector volumes. 
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Figure 4- The superposition of theoretical proton rate at solar minimum (dot- 
ted line) with the sum of five simulations with GSPM (continuous line) . GSPM 
reproduces very well the GCR spectral energy distribution. 



6.1. Preliminary Tests 

To test the correctness of our results we first made some preliminary 
tests. 

Tlie goodness of tlie internal randomness of GSPM was checked as 
follows. In Figure 4 both the theoretical (CREME94) and the G4 output 
of the GSPM energy distributions are shown. Five simulations with 
5x10^ particles each were run and the composite output is compared 
to the theoretical normalized one (see Figure 3). 

6.2. The r^sec/f^prim RATIO 

Table IV reports the number of secondary particles produced on average 
per primary event. The ^^sec/?^prim ratio is clearly increasing with the 
atomic number of the primary particle and the high and low stressed 
PhC have rather the same dependence. Unfortunately such a ratio is 
strictly dependent on the chosen physics list and on the reliability of the 
G4 models. G4 has been tested against nuclear data for nuclei lighter 
than or equal to C, therefore the ratios results on N and O must be 
taken with caution. 



6. Results from GCR simulations 
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Table IV. The ngec/^prim ratio from the simulation of GCRs radiation 

on the photoconductors. In the first column the GCR type; in column 
^2 and ^3 the average of secondary events generated by a single 
primary event for the high stressed PhC in case of solar min and 
max irradiation; in column #4 and #5 same as #2 and #3, but for 
the low stressed PhC. 



High Stressed Low Stressed 

Solar min Solar max Solar min Solar max 



GCR 


(Tisec ) / Tiprim 




(^isec ) /^Iprim 


{flsec} / Tiprim 


H 


3.302±0.122 


4.488±0.157 


3.272±0.088 


4.487±0.127 


He 


6.569±0.265 


10.612±0.279 


8.672±0.241 


10.385±0.335 


Li 


14.665±0.744 


17.409±0.416 


14.833±0.269 


17.394±0.493 


C 


35.859±0.818 


45.489±0.766 


36.473±0.913 


45.844±1.105 


N 


45.527±1.819 


57.420±0.913 


46.275±1.212 


57.658±1.015 


O 


55.863±1.394 


70.275± 1.188 


55.379±1.745 


72.563±1.234 



6.3. Deposited energy, definition of a glitch 

A glitch on the photoconductor is an unexpected voltage jump during 
the integration ramp. Along a single ramp, that lasts 300 ms, the 
voltage with a 64 readout sequences decreases monotonically. When 
the voltage jump exceeds by 4 - 6 times (in sigma units) the mean 
jump we have a glitch. 

The exact electronic behaviour is not reproducible with Geant4, 
we only collected the energy depositions on the pixels. The energy 
deposited on the pixels is transformed (Groenewegen 2004) in detector 
signal (Volt): 

AE{MeV) = C/W{Volt)Eg/{Rhv/ri) (1) 

where the photoconductive gain R = 30 A/W, at A=170 jim. (rj = 0.3, 
with a large error bar); Eg is the energy gap, and it is equal to 2.9 cV, 
and C is the detector capacity. For C = 3 pF, we get: Ay = 0.0134AE'. 

Impacts on pixels are identified and used to compute the deposited 
energy. Tables^ V and VI report the mean deposited energy (in V) for 
both (high and low) stressed PhC for shielded simulations, i.e. those 
simulating the effect of the spacecraft. 



® Tables report the sum of the voltage jumps normalized to a second of GCR 
impacting. This value must be corrected by the readout and ramp lasting. 
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Table V. Mean number of hits and energy deposition expressed in MeV, for 
the high stressed pixels with spacecraft shielding, secondary events at solar 
minimum and maximum. 







High Stressed pixel 




GCR 


Solar 


min 


Solar 


max 




in) 


{dEdep}{MeV) 


in) 


{dEdep}{MeV) 


H 


101.41±5.68 


60.21±3.19 


50.86±0.84 


27.35±1.83 


He 


17.46±0.66 


23.72±0.88 


8.77±0.23 


9.31±0.09 


Li 


0.019±0.001 


0.32±0.01 


0.089±0.004 


0.12±0.01 


C 


2.24±0.07 


4.64±0.18 


1.15±0.04 


1.86±0.09 


N 


0.76±0.02 


1.54±0.06 


0.39±0.01 


0.64±0.01 


O 


3.23±0.09 


7.07±0.27 


1.74±0.04 


2.89±0.15 



Table VI. Same as Table V for the low stressed pixels. 







Low Stressed pixel 




GCR 


Solar 


min 


Solar 


ma:K 




in) 


{dEdep){MeV) 


in) 


(dEdep){MeV) 


H 


105.29±5.63 


61.74±4.42 


47.98±0.95 


24.80±1.54 


He 


17.59±0.49 


24.21±0.99 


8.66±0.37 


9.41±0.59 


Li 


0.019±0.003 


0.332±0.001 


0.087±0.001 


0.119±0.003 


C 


2.22±0.05 


4.62±0.05 


1.18±0.03 


1.90±0.08 


N 


0.78±0.01 


1.63±0.09 


0.399±0.001 


0.64±0.03 


O 


3.33±0.05 


7.36±0.21 


1.72±0.04 


2.85±0.19 



6.3.1. High stressed Photoconductor 

We considered first the effects on the high stressed PhC of the primary 
particle impacts on the sliicld and the detector, and the production of 
secondary events from the shield and the detector. We ran five sim- 
ulations for each nucleus. Each simulation considers 2x10^ particles. 
Results (see Table V) have been normalized to 1 s of irradiation. 



6.3.2. Low Stressed Photoconductor 

We ran five simulations for each nucleus. Each simulation considers 
2x10^ particles. Results (see Table VI) have been normalized to 1 s of 
irradiation. 



gegapaper2_rev.tex; 5/02/2008; 12:01; p. 13 



14 



Bongardo et al. 



Tabic VII. First column: number of pixel affected by the negative 
glitch. Other columns: the mean percentage (5 simulations) of 
GCRs affecting a fixed pixel length, in the shield+secondary 
events configuration, for each GCRs type, and for the high 
stressed PhC. Similar values are found for low stressed PhC, 
and at solar maximum. 



Length 


P 


a 


(G4%) 
Li C 


N 





2 
3 
4 


q 0+2.7 
0.0_2 y 

U-O_o.6 




5.61^-^ 


7 c:+2.6 
'•^-2.6 


6.211-1 

1.811:1 

0.11°;? 


' -0-1.8 

1 Q+O.S 

1+°-^ 



6.4. Cross talks 

We checked whether in the G4 outputs there are coincidences {cross 
talks) between adjacent pixels. The number of pixels affected by one 
single GCR as resulting from G4 simulations are listed in Table VII for 
the high stressed PhC^. Values are given in percentage for each GCR 
type producing cross talks. The percent probability that nearby pixels 
are affected by a hit is lower than the values given in Table VII. 

6.5. Secondary Events 

The formation of secondary events depends strictly on the physics we 
considered. In case of proton GCRs we obtained a lot of secondary 
events, from photons and electrons to a large sample of particles (i.e. 
TT"*") to some nuclei (i.e. ^^Al). For the light nuclei, secondary events are 
basically photons or electrons. Ion inelastic scattering may also activate 
and generate pions, protons and electrons and smaller nuclei, but these 
events are much rarer. 

Figure 5 reports the number of impacts on the detector in a shielded 
and unshielded configuration (high stressed PhC). The shield does 
not act as the major source of secondary events and does not largely 
affect the pixels. Ge pixels are well embedded within other detector 
components masking the shield, the origin of the secondary. 

As expected, the nscc/''^prim ratio is high because of the large mass 
of the detector. The ratio is almost identical for high and low stressed 
PhC, since the two arrays are very similar. 



No significant dependence on solar activity or on stress status is found. 
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Figure 5. Cartesian projection of secondary events generation space for the 
high stressed PhC (from 5000 protons seed) for shielded (straight hne) and 
unshielded (dashed line) configurations. Upper panel; the projection along the 
Z axis, mid and lower panels projection along the Y and X-axis. 

7. Proton irradiation tests 

At the Light Ion Facility of the Centre de Recherches du Cyclotron 
of the University Catholique de Louvain La Neuve, Belgium, some 
tests were performed, intended to determine the glitch event rate on 
a single photoconductor module. The 2004 tests could not be repro- 
duced because of sudden unexpected responsivity changes which made 
impossible to follow the detector behaviour. In addition the external 
environment was not clearly defined and simulations could not usefully 
compared to measurements. Tests were repeated in April 2005 and the 
responsivity jumps were cured by heating the detectors (see Katterloher 
et al. 2005). Here we report on these latter tests. 

7.1. Test Specimen 

The test specimen foreseen for the proposed investigations under proton 
irradiation is the detector module FM#12 with a CRE of the Qualifi- 
cation Model type, mounted in the centre of the module. The detector 
is in the high stress configuration. The module stood in a liquid He 
dewar operating at a temperature of (1.85ib0.05) K. 
The dewar is made of three concentric cylinders of Al (called hereafter 
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Figure 6. The dewar, and a single stressed module, in the Louvain cyclotron 
configuration. In gray the three Al cylinders and the box around the module, 
in red the Cu cylinder foil. 

All, A12 and A13) and one made of Cu. All is 6 mm thick. A12 and 
A13 are 1.5 mm thick. The Cu cylinder is 0.5 mm thick. The module 
is placed inside a box of Al 4 mm thick [Al case) so that half the fore- 
optics stands outside of it (see Fig. 6). The module is not placed at the 
centre of the cylinders. 

The proton beam, before reaching the Ge:Ga crystals elements, pen- 
etrates two layers, one made of steel, 60 /xm thick, and one made of 
lead 0.12 mm thick. There is also a couple of steel collimators, that we 
modeled as a box, with a hole in its centre as large as the beam. 

During the test a primary circular (10 cm 0) proton beam of 70 MeV 
was used, with a beam length of 5.12 m between the diffusion foil and 
the Device Under Test (DUT). The beam reaches the dewar (and the 
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Figure 7. The UCL-CRC experimental hall. The detail of the DUT (here in 
the center of the plot) is in Figure 6. 

module) under an angle of 10° (see Figure 7). The DUT is placed into 
air. Due to air leakage, we set air cuts as large as we can, in order not 
to have a massive secondary event generation, that is 2.5 m. 

7.2. The G4 Simulations 

We first ran 5 simulations of 10000 protons each, in order to have an 
idea of what occurs to the beam once it crosses all the layers. We 
assume that the beam had a Gaussian error of 1 MeV. Before the first 
Al cylinder layer we had an energy of ~ 58.09 MeV, which disagrees 
with data by J. Cabrera, 63.74 MeV (Katterloher et al. 2005). This is 
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Tabic VIII. Summary of proton seed and photo-conduetor 
module setup during UCL-CRC experiments. Column # 1 
reports the identifier of the chosen setup, # 2 the value of 
the bias voltage of the FEE circuit. ^ 3 the duration of 
each integration ramp, # 4 the integrator capacitor, # 5 
the number of ramps, # 6 the number of repetition with the 
same measurement setup, # 7 the proton flux. 



Exp. 


BIAS 
(mV) 


tramps 

(ms) 


C 
(pF) 


N 

ramps 
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Flux 
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250 


0.43 
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10 




30 


250 


0.23 


512 
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10 


#111 


:!() 


2:){) 


0.2:! 


102 1 


18 


100 



due to the new G4 version we used (6.2.02^, gave ~63.70 MeV). Our 
modelled dewar is reliable and represents closely the dewar geometry. 
We do not have reasons to believe the Geant4 fails in this case and we 
set the energy hitting the module at 15.22±1.95 MeV (bottom right 
panel of Figure 8). In Figure 8 we plot the degradation of the proton 
energy along the path inside the dewar. 

What we would like to stress here is that 10000 protons originate 
roughly 22300 secondary events. This is mainly due to the air leak- 
age. After the Al box, ~610 protons survive. But, after the Al box, 
there are also saSSO secondary events. The spectral range and type of 
these secondary events is strictly dependent of the reliability of the G4 
models. 

In Table VIII we give the relevant information on the UCL-CRC 
experiment as performed at the Louvain cyclotron (Katterloher et al. 
2005). 



^ Furthermore the 6.7.01 version gave a very rare crash, probably due to a bug 
in the Hadronic Physics for the High Precision neutrons: with this new version we 
switched off such a physics. 
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Figure 8. Degradation of the proton energy across the dewar layers: from the 
top-left panel, the energy of the incoming proton flux, the flux emerging from 
the polystyrene foil, after the Aluminum shields, after the Cu case and the 
final detector case. The last panel refers to the energy deposition on the Ge:Ga 
material. 



7.3. The Geant4 results 



We ran the same experiments under the G4 tool. Results are summa- 
rized in Table IX. In Fig. 9 we plot the energy deposition onto the 16 
pixels in the HI test. 

Here wc found 3 peaks. The major one, is clearly due the primary 
proton beaming. The one at low energies, is due to the secondary events 
generated along the tracking of the protons along all the components 
we described before. The peak at higher energy could be due to the fact 
that the pixel are inside a non uniform cavity: due to the inclination of 
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Figure 9. Energy deposition for the HI test. 

Table IX. The UCL-CRC experiment outputs with Geant4. 



Exp. 




Flux 


^rad 




(events) 


Rate 








(cm-^s-^) 


(s) 
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253.20 


3.55 
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10 


10 
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100480 


129.10 


3.62 


#H1 


18 


18 


400 
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8038400 


9948.06 


139.58 


the b 


earn there is a 


part of it that 


cross 


the cavity 


in its thinner part, 



that is: pixel are hit by more energetic protons. 
7.4. Comparison with measurements 

A preUminary analysis of these data was made by Groenewegen & 
Royer (2005) . Here we report the results in energy; the detector output 
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signal in Volts can be retrieved from equation 1. In the report by 
Groenewegen h Royer (2005) the conversion value between MeV and 
Volt is taken to be 1.34. 

7.4.1. Events and Rates 

The number of events observed in files L26-L29 were 755, corresponding 
to a rate of 3.1 s~^cm~'^. If corrected for the efficiency of the detec- 
tion algorithm it becomes 3.9 s~^cm~'^. For such a dataset, G4 finds 
853 events and a rate of 3.53. The number of events observed in files 
II3-H6 is 29008, 04 has 29430 events. Within the uncertainties these 
values are similar. The fundamental physical processes occurring in this 
experiment axe therefore reproduced. 

7.4.2. Deposited Energy 

In Figures 9 and 10 we plot the distribution of the predicted deposited 
energy for the files L26, L27, L28 and L29 and those at large proton 
flux HI in MeV. The energy distribution of the glitches, i.e. that of the 
number of hits reaching the detector module, shows that many glitches 
occur at low energy and the distribution is asymmetric with a peak 
at Ri7 MeV. The distributions are bimodal and have a second peak 
around 8 MeV (sec §7.3) The energy glitch distribution is, as expected, 
the same in both experiments (HI and L26-L29) with low proton and 
high proton fluxes (Figures 9 and 10). What is changing is only the 
number of hits which in the second case is hundreds times lower. In 
fact, the proton energy is the same, in both configurations, only the 
beam intensity is changing. 

Figure 11 shows, in the upper panel, the measured glitch distribu- 
tion in Volt, and in the lower panel the high energy tail of the glitch 
distribution shown in Figure 10. We apply here the same deglitching 
algorithm and get rid of the low energy tail of the hits. This figure, 
i.e. the glitch detection, depends strictly on the algorithm used to de- 
glitch the data (Table 9 of Groenewegen &: Royer 2005). If, for instance, 
we applied a cutoff value around Ri6 MeV the distribution of glitches 
would be more similar to the observed one. This latter has a peak at 
low energy, with a rather long tail towards high energies (Groenewegen 
& Royer 2005). The simulated distribution gave first the tail, then the 
two peaks and is much narrower (cfr. §7.3). 

To compare the two distributions a conversion factor to the x- 
axis values (either to MeV or to Volt) must be applied. This can 
be accomplished using eq. (10). The resulting value is, in our case, 
AE{MeV) = 6.71AF(F). This value was computed from the measure- 
ments taking into account the responsivity changes under irradiating 
conditions (Lothar Barl, private communication). The two distributions 
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Figure 10. The energy distribution of glitches for files L26-L29, with Geant4. 

however differ substantially, the simulated one being much narrower. 
The difference could be due to a large number of secondary events 
produced during the tests by some unknown components we are not 
aware of. This component should have a large energy and it is unlikely 
that it occurs but this hypothesis cannot be discarded. Another pos- 
sible explanation is that the conversion factor from MeV to Volt is a 
function of the energy. Although plausible, at present we cannot prove 
this possibility. 



7.4.3. Boundary Effect 

Figure 12 shows the number of detected glitches per detector: a clear 
boundary effect (i.e. external pixels of the module were under-hit with 
respect to the central ones) is seen, and is due to a differential incident 
proton flux with respect to the detector position (Figure 12 in Kat- 
terloher et al. 2005). We find a difference of the hit numbers of 20% 
while the measured beam intensity difference between the central pixel 
and the outer ones is 10%. This is a geometrical effect and does not 
correspond to a different behaviour of the single chip. 



gegapaper2_rev.tex; 5/02/2008; 12:01; p. 22 



GCRs impact on PACS photoconductor with Geant4 



23 




Figure 11. Upper panel; the observed distribution of glitches once a deglitching 
algorithm has been applied. Lower panel: the G4 output of the same measuring 
setup and the same deglitching algorithm. 
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Figure 12. The boundary effect. Number of G4 simulated glitches per detector 
(pixel) in the HI setup. 
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8. Discussion 

8.1. The photoconductors in the space environment 

Our simulations of the effects of Galactic Cosmic Ray impacts on the 
photoconductor arrays show that both the satellite plus the cryostat 
(the shield) and the detector act as source of secondary events, affecting 
the detector performance. Secondary event rates originated within the 
detector and from the shield are of comparable intensity. The impacts 
deposit energy on each photoconductor pixel but do not affect the be- 
haviour of nearby pixels. These latter are hit with a probability always 
lower than 7%. 

The energy deposited produces a spike which can be hundreds times 
larger than the noise. The present simulations are not able to follow 
the temporal behaviour of each pixel and cannot be used to determine 
the shape of the output signal function of time. 

8.2. The photoconductors in the lab test 

We have simulated the experiment carried out at UCL-CRC and com- 
pared the simulation outputs to the measurements reported in Groe- 
newegen & Royer (2005). We find similar rates and events (see Tables 
X) on each pixel. The simulated energy distribution differs substantially 
from that of the measurements. 

We tried to ascribe these differences to the following causes: 

1. As far as the event number and rate is concerned, it may be possible 
that the chosen cuts are rather large. But due to the complexity of 
the experimental hall it is hard to tune them optimally. We would 
need information on the intermediate passage of the beam through 
matter. A simple correction could be a slight increase of the Ge cuts. 
A clear benchmark would be the comparison between the energy 
distributions. 

2. The presence of a primary and a secondary peak seems to be the 
only feature common between experimental and simulated data. 
What it is not clear is why, due to the crossing through matter, the 
Gaussian beam is (should be) transformed as in Figure 8, that is 
first a tail, then a peak, whereas experimental data favour the tail 
next to the secondary peak. Either Ge has some physical properties 
which are not included in the G4 Physics List, or the deglitching 
algorithm used is not correct. 
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Table X. Measured versus predicted values. 





(Groenewegen & Royer 2005) 


Geant4 


Files 


Events 


Events 


L26-L29 


755 


853 


H3-H6 


29008 


29430 
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